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Quenched Light Hadron Spectroscopy: Comparing the Wilson and 
0(a) -Improved Fermion Actions 
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The King's Buildings, Mayfield Road, Edinburgh, Scotland 

We have studied the light hadron spectrum and decay constants for quenched QCD at /3 = 6.2 on a 24^ X 48 
lattice. We compare the results obtained using a nearest-neighbour 0(a)-improved ("clover") fermion action with 
those obtained using the standard Wilson fermion action on the same gauge configurations. For pseudoscalar 
meson masses in the range 330-800 MeV, we find no significant difference between the results for the two actions. 
The scales obtained from the string tension and mesonic sector are consistent, but higher than those derived from 
baryon masses. The ratio of the pseudoscalar decay constant to the vector meson mass is roughly independent of 
quark mass as observed experimentally. 



1. INTRODUCTION 

1.1. O(a)-Improved Fermion Action 

Amongst the many systematic errors involved 
in the numerical simulation of QCD are those due 
to the non-zero lattice spacing, a. The standard 
Wilson pure gauge action differs from the contin- 
uum action by terms of 0{a^), while the Wilson 
fermion action introduces a term of 0{a) in order 
to avoid the fermion doubling problem. Sheik- 
holeslami and Wohlert proposed a form of the 
fermion action which was shown by Heatlie et 
al. to remove terms of 0{a) from matrix el- 
ements in perturbation theory. The improvement 
procedure described in ||^ involves rotating the 
fermion fields and adding a local "clover" ||] term 
to the Wilson action, i.e., 
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q{x)F^,^{x)a^uq{x), (1) 



where F^i, is a lattice definition of the field 
strength tensor. 

In this paper, we present results obtained for 
the light hadron spectrum on a 24? x 48 lattice 
at /3 = 6.2 using both the standard Wilson action 
and this 0(a)-improved clover action. 



1.2. Computational Details 

The simulations were performed on the Meiko 
i860 Computing Surface sited in Edinburgh. This 
consists of 64 nodes and has a sustained perfor- 
mance for QCD simulations of 1 to 1.5 Gigaflops. 
The gauge configurations were generated using 
Hybrid Overrelaxation Q. The overrelaxation 
was performed on a sequence of 3 SU{2) sub- 
groups of SU{3). As this is a micro-canonical up- 
date, we follow every 5 overrelaxed sweeps with a 
Cabibbo-Marinari heatbath ||] sweep to preserve 
ergodicity. We allow 16800 sweeps for equilibra- 
tion and separate configurations by 2400 sweeps. 

The comparison of the results for the two 
fermion actions is performed on the same 18 con- 
figurations. For each action, we have propaga- 
tors at 5 different values of the hopping parame- 
ter, K, chosen to give roughly similar pseudoscalar 
masses in the range 330 — > 800 MeV. As reported 
in 1^ , we found that the correlators obtained us- 
ing the clover action were somewhat noisier and 
so we have smeared the clover propagators at the 
sink to improve the signal. We also present re- 
sults for the clover action on 36 configurations 
with 3 K values and a local sink. 

The most computationally expensive part of 
this study was calculating the quark propagators. 
We find that red-black preconditioning improves 
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the efHciency by a factor of 3 and that, for the 
range of pseudoscalar masses investigated, Mini- 
mal Residual is typically twice as efficient as Con- 
jugate Gradient. In agreement with the study re- 
ported by Hockney in ^ , we see only a small gain 
from over-relaxation for MR at /3 = 6.2. 

2. RESULTS 

2.1. Fitting Procedure 

Limited statistics and the amount of supercom- 
puter time used imply that it is essential to per- 
form a careful analysis of the data. All our prop- 
agator fits take into account correlations between 
timeslices. We require that the fitting region, 
[imin, imax], satisfics the Condition that small per- 
turbations of tmin and tmax givc the same mass 
within errors. We find that the intervals [12, 16] 
for local sinks and [9, 13] for smeared sinks are 
satisfactory. The correlated x^/dof varies be- 
tween 0.3 and 4, which indicates both that we 
are taking correlations into account and that we 
are getting reasonable fits. Unlike the QCD- 
PAX collaboration Q], we find that increasing 
^min does not significantly change the masses ob- 
tained. However, we are only able to increase imin 
slightly before our signal becomes noisy. 

The errors are estimated using a bootstrap pro- 
cedure. iVcfgs configurations are chosen from the 
ensemble of A'cfgs configurations, allowing repe- 
tition. The same analysis is performed on this 
bootstrap sample as on the original A'cfgs config- 
urations. We repeat this typically 1000 times and 
determine the errors such that the middle 68% of 
our bootstrap samples lie within the error bars. 

2.2. Hadron Masses 

Throughout the range of quark mass used, we 
find that mj^g is linear with 1/2k, in accordance 
with PCAC. From correlated fits on 18 configu- 
rations, we obtain the following values of Kcrit- 
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The Edinburgh plot for the data is shown in 
Fig. |l]. The results for the two actions are broadly 
consistent, showing a trend towards the physical 
value of mjv/my with decreasing mps- 
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Figure 1. Edinburgh plot for both actions. The 
physical and heavy quark points are marked by 
the left-hand and right-hand □ respectively. 



For the comparison of the two actions, we have 
concentrated on quantities which might be more 
sensitive to 0{a) effects, such as mass splittings. 
In Fig. H, we plot the difference between the 
squares of the vector and pseudoscalar masses as 
a function of the square of the pseudoscalar mass 
in lattice units. Experimentally, this quantity is 
almost independent of the quark masses for both 
heavy-light and light-light mesons. We denote 
the range of experimental values by the horizon- 
tal lines, setting the scale using the string tension 
{i.e., a'^ = 2.73 tl GeV |]). In the chiral limit, 
the splitting for the Wilson action is compara- 
ble to the experimental values, while, at heavier 
quark mass, it is considerably smaller. For the 
clover action, there is a much weaker dependence 
on the quark mass and the values are roughly 
consistent with experiment. 

The differences between the A and nucleon 
masses are shown in Fig. ^. The leftmost points 
are obtained from the chiral extrapolations of the 
individual masses. These values, as well as those 
from lighter quark masses are in good agreement 
with the experimental value shown. 
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The experimental range is 



indicated by the horizontal lines. 



Figure 3. Baryon mass splitting. The experimen- 
tal value is marked with a □. 



In Table 0, we present various determinations 
of the lattice scale. The values obtained from 
hadron masses are from correlated linear extrap- 
olations to the chiral limit. The scales obtained 
from TUp agree well with the scale from the string 
tension, whereas those from the baryon masses, 
while self-consistent, are lower. 



Table 1 

Scales determined from various quantities. 

GeV 



Action 

^cfgs 



Wilson 
18 



clover 
18 36 



rriN 
mA 



2.77 

2.3911° 



-15 



2.57 
2.45 
2.33 



+16 
-16 
+18 
-11 
+10 
-10 



2.68 til 

+10 
-15 
+14 
-13 



2.38 
2.36 



string tension 



2.73 



+ 5 

- 5 



2.3. Decay Constants 

The results presented for the decay constants 
below are all obtained using propagators with lo- 
cal sinks. 

2.3.1. Pseudoscalar Decay Constant 

This is defined through the matrix element of 
the axial current 

Za{0\ Ai \PS) 



fps = ZaJps = 



mps 



(3) 



and our normalisation is such that the physical 
value is 132 MeV. We calculate fpg using 



fps 



1 



mps 



{PP) 



(4) 



with = qjfj.j5q and P = qj^q. We fix mps 
and (0| P \PS) from the pseudoscalar propaga- 
tor. In perturbation theory, the lattice renor- 
malisation factors are — 1 — QA3g^ and 

Z^ ~ 1 — 0. 02(?^ Using the bare coupling gives 
Zj ~ 0.87 and Z^ ~ 0.98, whereas, if we use 
the "effective coupling" proposed by Lepage and 
Mackenzie §, we get Zj ~ 0.78 and Z'x ~ 0.97. 

We have investigated the dimensionless quan- 
tity fps/'m-V: and find it varies only slowly with 
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quark mass, in agreement with the experimental 
observation that j-^jmp (0.17) is approximately 
the same as JkI'^.k* (0.18). In Table |[ we show 
the values obtained from linear extrapolations to 
the chiral limit. For both actions, the chiral ex- 
trapolations are roughly consistent with experi- 
ment. The fact that the data from the clover ac- 
tion has larger statistical errors is compensated 
by the smaller uncertainty in the lattice renor- 
malisation factor, . This means that we can, 
in principle, obtain more accurate predictions us- 
ing the clover action by increasing statistics. 



Table 3 

l//y as a function of n. 



Wilson 



clover 



!//• 



1// 



0.15100 


0.40 


+ 1 

-1 


0.14144 


0.33 


+ 1 

— 2 


0.15200 


0.43 


+ 1 

-1 


0.14226 


0.36 


+ 2 
-1 


0.15230 


0.44 


+ 1 

-1 


0.14244 


0.36 


+ 2 
-1 


0.15260 


0.45 


+ 1 

-1 


0.14262 


0.37 


+ 2 
-1 


0.15290 


0.47 


+ 2 
-1 


0.14280 


0.36 


+ 2 
-1 


'^crit 


0.47 


+ 2 
-1 


'^crit 


0.38 


-1-3 
-1 



Table 2 



Chiral extrapolations of fps/fny- 



Action 


Wilson 


clover 




18 


18 36 


fps/^v 


0.21 


0.13 0.15 




0.18^3 


0.12^2 0.15^2 




0.16^3 


0.12+^ 0.15^2 



2.3.2. Vector Decay Constant 

We define fv through 

i^^-'ir^^-^mW) (5) 

Jv Jv "'y 

with Vi — q^iQ- Table ^ shows the lattice results 
for both actions from 18 configurations. The val- 
ues depend only weakly on quark mass for both 
actions. The results obtained from a linear ex- 
trapolation to the chiral limit are given in the fi- 
nal row. Using the bare ("effective") coupling, 
these correspond to 1//^ = 0.39 +i (0.33 1^) 
and = 0.34+^ (0.31 tl)- These are similar 

to the experimental value, 1/ fp = 0.28(1). Using 
the conserved vector current would remove the 
uncertainty in the normalisation for the Wilson 
action. 



3. CONCLUSIONS 

For both actions, mass splittings and decay 
constants typically agree with experiment within 
errors. We find broadly consistent scales from dif- 
ferent quantities, though the baryon scales appear 
somewhat low. 

We have shown that the clover action can be 
used to study light hadrons. The computational 
overhead is small, except that the increased sta- 
tistical fluctuations mean that we require smear- 
ing or more configurations for comparable er- 
rors. We see a few hints of improvement, though, 
in general, the results are similar to those from 
the standard Wilson action. This suggests that, 
at P — 6.2, 0{a) effects may be small for light 
hadrons. 
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